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Abstract 

The existence of macroscopic regions with antibaryon excess in the baryon asymmetric 
Universe with general baryon excess is the possible consequence of practically all models 
of baryosynthesis. Diffusion of matter and antimatter to the border of antimatter domains 
defines the minimal scale of the antimatter domains surviving to the present time. A model 
of diffused antiworld is considered, in which the density within the surviving antimatter 
domains is too low to form gravitationally bound objects. The possibility to test this model 
by measurements of cosmic gamma ray fluxes is discussed. The expected gamma ray flux is 
found to be acceptable for modern cosmic gamma ray detectors and for those planned for 
the near future. 

It was recently stated HJthat the existence of antimatter domains with scale up to 1000 Mpc in 
a baryon-symmetric Universe cannot escape contradiction with the observed gamma background. 
This argument does not exclude the case when the Universe is globally baryon asymmetric and 
the antibaryons within antimatter domains contribute a small fraction of the total baryon charge. 
On the other hand, it was shown that the existence of antimatter domains in the baryon 
dominated Universe is a profound signature for the origin and evolution of baryon matter and its 
inhomogeneity. Depending on its parameters, the mechanism of inhomogeneous baryosynthesis 
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can lead to both high- and low-density antibaryon domains. According to 0, high-density domains 
can evolve into antimatter stellar objects, so that a globular cluster of antimatter stars can exist 
in our Galaxy, what may be tested in the cosmic searches for antimatter planned for the near 
future. On the other hand, as it was mentioned in [|3|, low-density antibaryon domains cannot 
evolve into gravitationally bound objects. The case of such "diffused antiworld" is considered in 
the present paper. 

A simple modification of the baryogenesis scenario will lead to formation of domains with 
different signs of baryon asymmetry. To achieving this, the condition should be satisfied that 
C- and CP- violation must have different signs in different space points . For instance, this can 
occur if there exist two sources of CP violation, a spontaneous and an explicit one. Spontaneous 
symmetry breaking naturally proceeded through the second order phase transition which started 
after the inflationary stage ||. One would expect that CP- odd amplitudes are not universal, space- 
independent quantities, but have different signs and values in different space points: (fi s = 4> s (x). 
If the amplitude of spontaneous symmetry breaking is small in comparison with the explicit one, 
<p s < <p e , there would be relatively small fluctuations in the baryonic number density around 
the uniform baryonic background. In the opposite case, <p e < <p s , the amplitude of fluctuations 
is large with respect to the baryonic background density. One can see that if the domains with 
high antibaryonic density are formed in the second-order phase transition, their elementary size 
li at the moment of formation is determined by U ~ 1/(AT C ) ||, where T c is the Ginzburg 
critical temperature at which the phase transition takes place and A is the self-interaction coupling 
constant of the field which breaks CP symmetry. In this case different elementary domains would 
expand together with the Universe and now their size would reach the value Iq ~ li(T c /To) = 
1/(AT ) ~ lCT 5 cm/A, where T is the present temperature of the background radiation. 

Another possible mechanism for generation of large fluctuations of baryon asymmetry in space 
is based on the model of baryogenesis with baryonic charged condensate. At first, the baryon 
asymmetry is accumulated in the form of a condensate of the scalar superpartner of the colorless 
and electrically neutral combination of quark and lepton fields (X) [|5|. This condensate might 
be formed during the inflationary stage if baryonic and leptonic charges were not conserved and 
the potential U(X) had a flat direction. The subsequent decay of this condensate would result 
in a considerable baryon asymmetry. Depending upon the initial conditions in the different space 
points, it is possible to produce baryons or antibaryons in the decay of the condensate (X) 0. The 
characteristic size of an elementary domain with a definite sign of baryonic charge was estimated 
in 0. At the end of inflation it is equal to li ~ H~ x exp(A x 1 ^ 2 ), where H is the Hubble parameter 
at the end of inflation and Ax is the constant of self-interaction of the field X. 

The next possibility to generate antimatter domains is connected to the existence of topological 
defects of the CP-violating axionic field. The existence of the vacuum angle 9 leads to strong CP 
nonconservation The famous solution of the strong CP problem is the introduction of the 
Peccei-Quinn (PQ) symmetry together with an axion field || [|. In this solution, 9 has the 
dynamical meaning of axion field amplitude, and goes to zero for zero vacuum expectation value 
of this amplitude. However, in the early Universe the effective vacuum angle 9 varies in the range 
[0, 27r], and this could imply a significant role of strong CP violation in baryogenesis. The point 
is that, in the standard model of invisible axion, one of the consequences of the effective vacuum 
angle 9 is the existence of P- and CP- odd condensates. Consider the cosmological properties of 
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the N = 1 axion model. Here N = 1 means that the axion model has a unique vacuum. The 
following potential embodies the qualitative features of the model of interest to us: 

L = \d^d^ - - A (<P + <P - F 2 a f + A 4 (l - costf) (1) 

where = (f>i + i<p2 is a complex scalar field. For A = this model has a U(l) global symmetry 
under which <j)(x) i— > exp(i9)<j)(x). This U(1)pq symmetry is spontaneously broken by the vacuum 
expectation value ((f)) = F a exp(i9). The associated Nambu-Goldstone boson is the axion. The 
last term in Eq.(l) represents the non-perturbative QCD effects that give the axion its mass 
m a . The model has global string solutions for A = 0. A straight global string along the z- axis 
is the static configuration (0) = F a f(p)exp(i9), where (z,p,9) are cylindrical coordinates. The 
surprising result that emerges from numerical simulations is that most of the energy in the string 
network corresponds to infinite strings. Precisely, the density in infinite strings was estimated to 
be the 80% of the total string density [T(| . 



In the early Universe, the standard axion is essentially massless in the time interval from the PQ 
phase transition at temperature of order F a to the QCD phase transition at temperature ~ lGeV. 
During that epoch, axion strings are present as topological defects. At first the strings are frozen 
in the plasma stretched by the Hubble expansion. However the plasma becomes more dilute with 
time and, when the temperature decreases below T w 2 ■ 10 7 GeV(F a /10 12 GeV) 2 the strings 
move freely and decay efficiently to axions. Since this decay mechanism is very effective, there 
is approximately one string per horizon from the temperature T to the temperature T\ ~ lGeV, 
when axion acquires mass. The existence of linear topological defects implies that the effective 
angle 9 varies by 2ir around the string, and this fact might be responsible for spatial modulation 
of the baryonic charge density in the Universe. Let us suppose, according to ||, that, owing to 
the spatial variation of the CP -violating axionic phase, the baryon excess is also space dependent 
and given by 

B(x) = A + bsin9(x). (2) 

Here A is the baryon excess induced by explicit CP-violation processes, which gives rise to global 
baryon asymmetry of the Universe, and b is the measure of the baryon asymmetry dynamically 
induced by axion. The antibaryon excess is generated within a wedge, which is symmetrical 
relative to a sheet along the direction 9 = 3tt/2, and the facet of this wedge is connected to the 
axionic string. The intersection of such wedges gives a domain with antibaryonic excess. Thus the 
typical size of an elementary domain will be as large as the horizon scale Zj ~ mpi/T^ s , where mpi 
is the Plank mass and Tps is the temperature of baryogenesis. Note that the defect formation can 



be triggered by inflation JT2[. If it takes place during the end of the inflation epoch, the defects 
will not be completely diluted away but can increase the typical distance between strings and, 
consequently, the typical size of an elementary domain. 

In any case all the above-mentioned mechanisms are able to generate domains of antimatter of 
any baryon density and on any scale. Such domains may consist of many elementary domains of 
small scale, as in the case of the standard phase transition, or it can appear like one moderately 
inflated elementary domain as in the case of the phase transition triggered by inflation. 

Let us consider the case of low-density antimatter regions, corresponding to the model of 
diffused antiworld. 
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When the density of antimatter pi within a domain is 3 orders of magnitude less than the 
baryon density pb (which we assume in the further discussion corresponding to flbh 2 = 0.1 , where 
h is the Hubble parameter normalized to lOOkms^ 1 Mpc" 1 ) , the cosmological nucleosynthesis in 
the period t ~ 1 — 10s results in a nontrivial chemical composition. For 10~ 4 < pi/pb < 10~ 3 
antideuterium is the dominant product. For smaller densities of antibaryons within domains 
no antinuclei are formed. At the density pi/ Pb < 10~ 4 , owing to the low antimatter density 
inside the domain, no recombination takes place at z ~ 1500 and the antimatter domain remains 
ionized after recombination in regions of baryonic matter. The radiation pressure and the energy 
density dominance within the domain suppress then the development of gravitational instability, 
so that the antimatter domains of sufficiently large size should be now clouds of ionized positron 
- antiproton plasma presumably situated in voids. 

Let us give a quantitative estimation of the surviving size and the observational effects of 
domains of a diffused antiworld. Let us assume that the Universe contains regions of very small 
antibaryon excess density. At temperatures above several MeV these regions cannot be strongly 
affected by the diffusion of surrounding particles, because their mean free path length is small 
enough. Therefore we will consider the evolution of these antibaryon domains at temperatures 
4 • 10 3 K < T < 10 9 K, when the Universe is radiation - dominated and contains mainly electrons, 
protons, photons and neutrinos. If the size of the antibaryon region is much greater than the 
mean free path of surrounding particles we can solve a one - dimensional problem assuming that 
the "initial" baryon density at T = 10 9 i^ is given by 

n fe (R,to) = | 0j x>Q (3) 
In this case the diffusion equation for baryons is 

-w = D{t) ^- anb (4) 

where D(t) is a diffusion coefficient. The last term in Eq.(4) takes into account the expansion 
of the Universe. Note that, since the antibaryon component is very small, we neglected it in the 
diffusion equation. 

Let us introduce a new variable r defined as the baryon to photon ratio, r = n 6 /n r Since the 
evolution of the photon density is given by the equation dn-Jdt = —an 7 , we can rewrite Eq.(4) 
in terms of r as 

< 6 > 

where 

r(R,t n ) = 




We note that in the problem under consideration photons are distributed uniformly in baryon 
domains as well as in the antibaryon ones. Due to electroneutrality of the media, the motion of 
protons is strongly related with the motion of electrons. The last ones interact with photons, and 
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this interaction defines the diffusion coefficient D(t) for protons. According to |13| this coefficient 
(in units where the Boltzmann constant kf, = 1) is given by 



D (t) « « 0.61 ■ 10 32 Z- 3 cm 2 / S , 



(6) 



where T 7 and p 7 are the temperature and the energy density of the radiation, respectively, c is 
the velocity of light, ot is the Thomson cross section, Z is the red shift which is related with the 
time t at this stage by t ~ 2.6 • 10 19 s/Z 2 . It is convenient to extract the t - dependence and to 
write D(t) in the following form 

D(t) = D(to) (£) 3/2 • (7) 

At the initial time to, corresponding to the temperature To = 10 9 K, D(to) is given according to 
Eq.(6) by 

D(t ) w 1.24- 10 6 cm 2 /s. (8) 



Solving Eq.(5) with the diffusion coefficient (6) we find 



r(R,t) 

where $ is the error function. 
Assuming t » t we get 
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or, in terms of temperature, 
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If we choose T = 4000K at the end of the radiation - dominated epoch, then 
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It follows from Eq.(12) that the average displacement of the boundary between baryon and an- 
tibaryon domains is about Ax ~ 0.2pc. Therefore, the primordial antibaryon regions of low 
density, which grow up to lpc or more at the end of the radiation - dominated stage, remain 
practically unaffected by the diffusion of the ordinary matter. 

Note that according to Eq.(9) in the limit t » to the motion of the boundary is determined by 
the final time t and does not depend on to- This is due to the fact that at the end of the radiation 
- dominated stage the particle number density drops significantly and as a consequence the mean 
free-path length increases, thus giving an important contribution to the boundary motion. 
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Below 4000K atoms are formed in baryon domains. Since the antibaryon density is assumed 
to be small enough (pi/pb < 10~ 4 ) in our approach, we can consider the flow of hydrogen atoms 
into antibaryon regions as a motion of free - streaming atoms. The physical distance traveled by 
the atoms after recombination until the present time t p is given by 

f*v v(t)dt , . 

d^a p 13 
P Jt rec a{t) y 1 

where v (t) is the average velocity of atoms and a(t) is the scale factor of the Universe. It is 
convenient to set a = ya p , where y = 1/(1 + Z) = = T p /T , and to rewrite Eq.(13) as 

d _ v(a)da = /-i v(y)dy 
Ja rec &<i Jarec/cip yy 

At the instant of recombination the velocity of atoms is given by the thermal value v rec ~ 



cJT rec /m, where m is the mass of the atom. After 4000^, the typical velocity of atoms v ~ p/m 



is red - shifted down as l/a ~ T. Therefore below 40007^ we have v ~ c(T p /m)Jm/T(l/y). 
Taking the equation 

{yf = y g p^ 2 ( 15 ) 

(where G is Newton's constant) and substituting for the density the expression p(y) = po/y 3 , in 
the case of a matter - dominated Universe we find 



d ~ / 3 CTp C — = / 3 j ™ (16) 

V 87rGp rn V T rec JT P /T rec y 3 / 2 y 8nGp y 

Tree 

Substituting p = p c , where p c = 1.88 • 10~ 29 h 2 = g/cm 3 is the critical density of the Universe, 
we obtain d ~ 3/h kpc. In the wide range of the Hubble constant h between 0.4 and 1 the free - 
streaming length of atoms will be of the order of several kpc. Therefore antibaryon regions of the 
same size will be filled at present by hydrogen atoms. 

A key observation to test the model of diffused antiworld could be the search for gamma rays 
from a boundary annihilation of antimatter and hydrogen atoms. 

Let us consider first the possibility of annihilation in antimatter domains filled with hydrogen 
atoms (the annihilation of matter - antimatter domains in baryon symmetric Universe was con- 
sidered in [0]). The equation for the number density of antiprotons which takes into account both 
the annihilation and the expansion of the Universe is given by 

^ = -{av)n b n- b - an- b . (17) 

In the limit n& >> ni we can neglect the variation of n b due to annihilation. Then , introducing 
r = rib/n-y, f = n&/n 7 , and solving Eq.(17), we find that at present time in antimatter domains 
filled with hydrogen atoms 

^p = ^re C exp[— / (av)n^rdt}. (18) 
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Since according to 
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at energies below lOeV the cross section of pH annihilation is given by 



((tv) « 2.7- 10- 9 cm 3 /s (19) 

we find that the integral in Eq.(18) is much greater than unity. So, no gamma radiation occurs 
at present from such regions because practically all antiprotons have been already annihilated. 

Therefore the radiation comes only from the narrow region of the boundary between matter 
and antimatter domains, where antibaryons have not annihilated yet. The width of this region is 
d ~ v At, where At is defined from the condition 

(av)n b dt « (av)n b At ~ 1, (20) 

-At 

and the velocity v is given by 



^ / ^ m (21) 



Substituting all numerical values, we find for the width of the region : 



T v T v 1 fl0- 7 cm- 3 \ . , 

d « cJ3-^^-— « 0.86 pc. (22) 

V T rec m (av)n b \ n b ) 



The gamma flux near Earth in this case is given by 



d$ dndN^ V 



dudVt dt duj 47rr^' 



(23) 



where dn/dt = {av)n b n b is the rate of annihilation per unit volume and per unit time; dN^/duo is 
the differential cross section for an inclusive gamma production; V = 4nR 2 d is the volume of the 
annihilating boundary of the diffused antiworld at present; R is the size of the diffused antiworld 
region; ta is the distance between Earth and the diffused antiworld region. 
Integrating Eq.(23) over photon energy we obtain 

where (iV 7 ) is the average number of photons per one act of annihilation (typically (iV 7 ) ~ 4) , 
(3 = n b /n b « 1. Finally, taking into account Eq.(20) we get 



The absence of absorption bands in the spectra of long-distance quasars indicates that the 
temperature of the intergalactic gas is high enough and that this gas is ionized. If this is the 
case we have to change the cross section in Eq.(19) by the corresponding cross section for pp 
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(av) w 6.5 • VQr X7 {v/c)cm A /s. (26) 
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However, for pp annihilation the integral in Eq.(18) will be much greater than unity, similarly 
to the pH case. Therefore expression (25) for the photon flux remains valid for pp annihilation 
because this flux does not depend explicitly on (cry). 

Up to now, no evidence for a major anisotropy in the gamma background has been observed. 
Below 1 GeV at high latitudes the diffused photon background is given by jnj 

« lO-WW- 1 . (27) 

According to a diffused antiworld can exist if /3 < 10 -4 . Therefore it follows from Eq.(25) 
that the possible gamma flux from the annihilation at the boundary of two worlds should be less 
than 10 -8 cm~ 2 s _1 sr _1 . This means that such diffused antiworlds could exist , in particular not 
far from our Galaxy, successfully avoiding to be detected. However, the development of modern 
detectors of gamma rays, like EGRET and AMS with the flux sensitivity up to 10~ 8 cm~ 2 s _1 sr _:L , 
gives a hope to detect diffused antiworlds. 

In this article we have found that the minimal scale of an antibaryon domain, which was not 
destroyed by annihilation, must be larger then 1 kpc at present. Independently of the conditions 
of creation of such an antiworld domain, the density within it will not increase due to gravitational 
instability, provided that two conditions are satisfied. Namely, this domain must be located in 
the void and the density of antimatter must be less than the critical density. It means that the 
space structure and geometrical form of the annihilation boundary of such a domain is conserved. 
In the case when the antiworld region of scale R consists of many antimatter domains of minimal 
annihilation scale d, the annihilation in such a region takes place at the border of each little domain. 
As a result the annihilation occurs in the whole volume of this region and, as a consequence, the 
gamma flux increases, owing both to the volume effect and the possibility for the antimatter 
density within a small domain to be as high as (3 m 1. This picture of volume annihilation in 
antiworld regions can be tested in searches for the angular variation of gamma flux in the angular 
range a ~ R/ta- 
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